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结合 纹理 复杂 度 和 JND 模型 的 图 像 水 印 算法 
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摘 要 : 针对 目前 基于 灰 度 共生 矩阵 相关 水 印 算法 存在 容量 较 小 的 问题 ， 提 出 一 种 将 纹理 复杂 度 和 DCT 域 JND 模型 
相 结合 的 图 像 水 印 算 法 。 首 先 将 原 图 像 分 成 子 块 ， 利 用 各 子 块 灰 度 共 生 短 阵 的 四 个 纹理 特性 计算 其 纹理 复杂 度 ， 并 据 
此 对 子 块 进行 排序 以 确定 水 印 瞪 入 的 子 块 位 置 后 对 原 图 像素 矩阵 进行 DCT 变换 ,结合 新 的 分 区 方式 计算 各 块 JND 值 ， 
根据 JND 值 以 及 新 的 谈 入 规则 确定 子 块 内 水 印 误 入 方式 。 算 法 有 效 考 虑 了 图 像 块 的 纹理 特性 和 人 眼 视觉 敏感 性 ， 加 强 
了 谈 入 水 印 后 的 图 像 质 量 ， 提 高 了 水 印 座 入 容量 。 实 验 结 果 表明 ， 所 提 方 法 在 谱 入 相同 容量 的 水 印 时 ， 图 像 平 均 峰 值 
信 噪 比 (PSNR) 较 现 有 方法 高 出 4.27% 左 右 ， 在 误 入 超过 原 方 法 容量 上 限 一 倍 的 水 印 时 ， 图 像 平 均 PSNR 仍 有 53.449 
8 dB。 
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Image watermarking algorithm combining texture complexity and JND model 


Li Shuzhi’, Long Xiangyu’, Deng Xiaohongb, Zhou Yongxina 
(a. School! of Information Engineering, b. School of Applied Science Jiangxi University of Science & Technology, Ganzhou 
Jiangxi 341000, China) 


Abstract: For a low capacity of the watermarking algorithm based on the gray scale symbiosis matrix, this paper would introduce 


the image watermark algorithm which combines texture complexity with DCT-domain JND model. Firstly, it divided the original 


的- image into several blocks, and the blocks the texture complexity would be calculated by four texture features of GLCM of each 


二 - block. The blocks would be sorted based on the texture complexity to locate the blocks of watermark embedding. After that, the 


prime matrix of original image would be transformed into DCT domain and, combined with new partition way, JND value of 
each block would be calculated. The mode of watermark embedding would be determined on the basis of JND value and new 
rules of embedding. The experimental results show that In the same capacity of watermark embedding, the average peak sinal- 
to-noise ratio(PSNR) of the image in this method is 4.27% higher than existing mathods, when embedding a watermark that is 
twice the upper limit of the capacity of the original method, the average PSNR of the image is still 53.4498 dB. 
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0 ”引言 等 人 四 提出 了 一 种 基于 DCT 变换 的 二 值 序列 分 组 数字 水 印 算 
法 , 利用 一 个 DCT 系数 租 入 多 比特 水 印信 息 , 减少 了 水 印 藤 入 

纹理 特征 分 析 最 初 应 用 于 遥感 图 像 分 析 03 ,由 于 纹理 特征 过程 中 需要 修改 的 系数 量 ， 提 高 了 水 印 算法 的 嵌入 容量 ， 实 现 

是 在 图 像 计算 中 经 过 量化 的 图 像 特 征 ， 通 过 对 其 分 析 能 够 得 到 ”了 水 印 嵌入 位 置 的 随机 选取 ， 并 对 典型 攻击 具有 重 棒 性 ，Shen 
图 像 有 效 的 宏观 和 微观 信息 ,所 以 纹理 特征 分 析 成 为 图 像 处 理 、 ”等 人 铝 提 出 了 一 种 基于 灰 度 共生 抢 阵 纹理 特征 的 无 源 拼 接 检 测 
数字 水 印 等 方面 的 重要 研究 途径 。 常 用 的 图 像 纹理 分 析 方法 有 ”方法 即 TF-GLCM。 在 TF-GLCM 中 , 对 不 同 的 块 进行 离散 余弦 


统计 分 析 、 结 构 分 析 、 模 型 分 析 和 频谱 分 析 。 而 统计 分 析 中 由 变换 ， 然 后 计算 GLCM， 以 充分 捕获 纹理 信息 与 图 像 像素 之 间 
HaralickD5 提出 的 灰 度 共生 和 矩阵 (gary level co-occurrence 的 空间 关系 ， 来 减少 特征 向 量 的 维 数 和 计算 复杂 度 。 李 淑 芝 等 
matrix,GLCM) 是 目前 统计 分 析 方 法 中 最 有 效 的 方法 之 一 。 张 浩 ”人 后 利用 GLCM 得 到 四 个 特征 参数 构建 了 一 个 计算 图 像 复杂 度 
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的 数学 模型 来 指导 水 印 的 圣 入 ， 该 方法 在 一 定 程度 上 提高 了 图 
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9 对 比 度 D。 对 比 度 可 以 用 来 描述 图 像 的 清晰 程度 ， 即 纹 


像 的 峰值 信 噪 比 (peak signal to noise ratio，PSNR),， 但 

嵌入 规则 的 局 限 性 而 限制 了 水 印 的 散 入 容量 ， 同 时 其 水 印 集 中 
在 

鲁 


芷 图 像 纹 理 最 复杂 的 128x128 子 块 上 骨 入 ， 影 响 了 图 像 质量 和 


与 此 同时 , 基于 DWT 或 DCT 域 的 JND 模型 中 也 逐渐 广泛 
的 应 用 于 数字 图 像 水 印 之 中 ， 李 智 等 人 [31 通过 考虑 人 眼 视 觉 掩 
蔽 特性 提出 了 一 种 DCT 和 DWT 的 综合 水 印 模型 , 最 小 可 觉 差 
(just noticeable difference，JND) 描述 人 眼 不 能 察觉 的 最 大 图 

像 失真 ， 考 虑 人 眼 能 够 接受 图 像 一 定 的 改变 而 不 被 察觉 。 在 图 
像 处 理 领 域 ，JND 能 用 来 检测 人 眼 对 图 像 中 不 同 区 域 失真 的 敏 
感性 。 近 年 来 ， 各 种 JND 模型 被 广泛 应 用 于 图 像 领域 。Hsu 等 
人 站 将 反 向 传播 神经 网 络 技术 和 JND 模型 结合 应 用 ， 水 印 目标 
系数 的 调整 受 JND 影响 和 BPNN 预测 ， 其 方案 能 承受 各 种 图 
像 处 理 攻击 ， 具 有 优异 的 鲁 棒 性 和 不 可 感知 性 ， 实 现 了 有 效 的 
盲 图 像 水 印 。 郑 明 鬼 等 人 09 为 了 提高 变换 域 JND 的 精度 , 在 计 
算 对 比 度 掩盖 因子 时 考虑 了 纹理 分 量 滤波 改善 了 ”JND 低估 问 
题 。Wu 等 人 0 在 计算 纹理 掩盖 时 更 深层 次 的 考虑 了 人 有 眼 对 图 
像 规则 区 域 和 非 规则 区 域 具有 不 同 的 敏感 性 ， 提 出 了 一 种 基于 
亮度 自 适应 与 结构 相似 性 的 JND 模型 。Hu 等 人 [9 为 了 增强 水 
印 的 不 可 感知 性 ， 利 用 JND 来 调节 参数 的 变化 幅度 ,提出 一 种 
基于 DCT 的 言 水 印 方案 , 将 二 进 制 信息 嵌入 图 像 ， 对 常见 攻击 
的 鲁 棒 性 具有 显 着 的 提高 。 叶 间 吕 提出 了 基于 人 类 视觉 特性 
HVS 的 DWT 域 盲 数字 水 印 算 法 , 同时 利用 JND 提高 了 数字 水 
印 的 不 可 感知 性 。 
综 上 所 述 ， 基 于 GLCM 的 纹理 分 析 和 变换 域 JND 模型 在 

水 印 的 嵌入 中 各 自 有 着 极为 广泛 的 应 用 前 景 ， 但 目前 还 没有 一 
种 方法 将 GLCM 与 DCT 域 JND 模型 结合 起 来 以 弥补 各 自 模 型 
的 不 足 ， 因 此 本 文 提 出 一 种 将 两 者 结合 的 水 印 岩 入 新 方法 。 本 
文 方法 优化 了 JND 模型 ,扩充 了 水 印 嵌 入 容量 ， 提 高 了 藤 入 水 
印 后 图 像 质量 。 


1 ”纹理 复杂 度 与 JND 衡量 
1.1 GLCM 复杂 度 计算 


克 度 共生 和 矩 阵 表示 的 是 图 像 中 像素 距离 和 角度 的 矩阵 函数 。 


GLCM 中 的 图 像 像素 的 距离 d=(d%dy)j 方 向 为 0 GE 
{0°%,45°,90°%135°)}, 即 d E{(0,d),(d.d),(d,0),(-d,d)?。 用 g(ij,q,0) 
表示 GLCM。 从 灰 度 共生 和 矩阵 中 提取 出 如 下 四 个 主要 的 特征 参 
数 : 


a) 能 量 J。 能 量 反映 了 图 像 分 布 均匀 程度 和 纹理 粗细 粒度 ， 


有 
Iag) (1) 
b) 业 。 博信 用 来 描述 图 像 所 包含 的 信 息 量 ， 
i j,d,0)og,g(i,j,d,0) O) 


i=l j=l 


时 的 清晰 度 ， 有 


D= > Yi-)) glijd,0) (3) 

d) 相关 性 COV。 相 关 性 可 以 用 来 描述 图 像 纹理 局 部 变化 
的 多 少 ， 有 

coOv=[ >( x j)g(i,j,d,0) -p41/ 00 (4) 


利用 均 方 误差 给 四 个 特征 参数 分 配 权 值 ， 可 以 得 出 计算 图 
像 纹 理 复 杂 度 四 的 计算 方式 如 下 : 


MSE, = (J,-— 7) (5) 
J= (J +J,+J;+],) (6) 
a = MSE 
DMSE, (7 
i=l 
J=0| xJI+0, xJ,+0, xJ,+0, xJ, (8) 
f=@ xJtw,xH+to,xD+to, xCOV (9) 


其 中 : 厂 、 有 J、J3、J4 分别 表示 GLCM 四 个 方向 上 的 能 量 ，O1、 
60.、63、04 是 根据 均 方 误差 分 配给 不 同方 向 上 同一 个 特征 参数 的 
权 值 。 根 据 式 (9)， 对 一 幅 512x512 标准 灰 度 图 像 Lena 按照 
32x32 分 块 并 计算 出 各 块 复杂 度 ， 实 验 结果 如 表 1 所 示 。 


表 1 部 分 图 像 的 图 像 复杂 度 信 息 


子 块 的 复杂 度 排名 ” 子 块 所 在 行 ” 子 块 所 在 列 ” 子 块 复杂 度 f 
1 193->224 ”449->480 54.34 
2 65->96 385->416 40.71 
3 129->160 481->512 32.72 
5 161->192 449->480 31.05 
10 225->256 417->448 24.12 
15 289->320 417->448 20.69 
25 1->32 1-32 18.25 
50 385->416 449->480 12.26 
在 得 到 各 块 复杂 度 排序 之 后 ， 将 据 此 结合 改进 的 DCT 域 


JND 模型 来 确定 嵌入 水 印 的 位 置 以 及 符 入 的 先后 顺序 。 

1.2 结合 纹理 特性 的 DCT 域 JND 值 的 计算 

基于 DCT 域 的 JND 模型 4 描述 为 空间 对 比 度 敏感 函数 、 

亮度 自 适应 因子 和 对 比 度 掩盖 因子 三 者 的 乘积 ， 关 系 式 如 下 : 
JND(i, D=J esp (is 7D) X A X Foalis )) (10) 

其 中 : Jcsr(ij) 为 空间 对 比 度 敏感 函数 ， 其 求解 公式 如 下 : 


5s、 (a 十 bo, Jexp(-c@,,) 
$9 r+(l—rcosp,) 


Jcsr(i,]) = (11) 
s 定义 为 集合 效应 ， 本 文 取 s=0.25, rf 取经 验 值 0.6，$,$ 
为 DCT 归 一 化 系数 ，9; 代表 相应 DCT 系数 的 方向 角 。 
Atum 为 亮度 自 适应 因子 ， 人 了 眼 对 图 像 亮 度 不 同 的 区 域 具 有 
不 同 亮度 自 适应 加 权 因 子 ，Ann 与 图 像 区 域 的 平均 亮度 工 存在 
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如 下 关系 : 按照 如 下 规则 来 划分 纹理 区 和 边缘 区 : 
一 - True， 纹理 区 

1+(60-1)/150,T < 60 oe (9) 
False， 边 缘 区 


A，14160 < I<170 (12) 


1+(T170)/425,T >170 


lum 


其 中 : 63、64 分 别 是 计算 图 像 复杂 度 中 对 比 度 和 相关 性 对 应 的 
权 值 。 根 据 上 述 构建 的 数学 模型 将 图 像 划分 成 平滑 区 、 纹 理 区 
Feontast(ij) 为 对 比 度 掩 盖 因 子 ， 考 虑 到 子 带 系数 自身 的 掩盖 ”和 边缘 区 。 根 据 图 像 区 域 的 分 类 , 不 同 区 域 具有 不 同 的 Feneas(ij) 


I 


效应 ， 最 终 得 到 如 下 关系 式 : 值 ,最 终 得 到 各 个 子 块 的 JND 值 .基于 GLCM 纹理 特性 的 DCT 
Fm (i,j) = 域 JND 模型 算法 如 算法 1 所 示 。 
w,(F 7) <16 且 处 于 平滑 区 或 边缘 区 算法 1 基于 GLCM 改进 的 DCT 域 JND 算法 
a (13) Input: 工 // 原 始 图 像 。 
yx wan "| rn) ,其 他 Output:M // 图 像 的 JND 值 矩 阵 。 


1: read(I); // 读 取 原 始 图 像 
在 图 像 不 同 的 纹理 区 域 Funms(ij) 中 y 的 值 不 同 ， 对 于 平 。 。 2:compute Jarfty 了) // 按 照 2.2 节 计算 公式 
滑 区 和 边缘 区 ， 人 眼 对 其 变化 较为 敏感 ， 所 以 式 中 y=1， 对 于 3:compute Am // 按 照 2.2 节 计 算 公式 


纹理 区 域 ， 由 于 人 有 眼 对 纹理 区 域 低频 系数 的 敏感 度 相 对 较 小 ， 4:num=512*512/8*8;// 将 图 像 分 成 num 个 8*8 块 
此 时 加 权 系 数 y=2.25， 而 高 频 系 数 的 加 权 w=1.25。 5:for each block B in {B,B,,…,B,, }// 对 每 个 子 块 操作 
对 于 一 幅 图 像 , 当 灰 度 共生 和 矩阵 中 元 素 集 中 分 布 时 , J 值 较 6:compute eo,; 


大 ， 表 明 图 像 是 比较 均匀 和 规则 的 变化 图 像 ， 反 之 则 较 小 。 妆 
灰 度 共生 矩阵 值 分 布 均匀 ， 也 即 图 像 近 于 随机 或 噪声 很 大 时 ， 


num num 
7: 已 =》e /num 》 a=>a /num ;3 
i=l i=l 


其 粹 值 较 小 ， 反 之 则 较 大 。 平 滑 区 是 一 种 比较 均匀 规则 的 纹理 8:if 。>E  ”// 区 域 类 型 判定 
模式 ， 所 含 信 息 量 较 非 平滑 区 少 ， 因 此 本 文通 过 焙 值 和 能 量 二 w=1 ;// 平 滑 区 
分 平滑 区 和 非 平滑 区 ， 具 体 过程 如 下 : else if 5>56 
首先 计算 出 整 副 图 像 的 平滑 度 E : wy =1;// 边 缘 区 
else if p22<16//h、1 分 别 为 子 块 中 各 位 置 的 行 、 列 
E=0, xJ+0, x| (14) a 
Ww =2.25 ;// 纹 理 区 低频 
其 中 : 6、6: 分 别 是 计算 图 像 复 杂 度 中 能 量 和 粹 对 应 的 权 值 ， else 
B Ww =1.25 ;纹理 区 高 频 
本 到 ; B 代表 图 像 分 成 的 块 数 ， 如 一 幅 512x512 的 图 像 ， end if 
B end for 
最 多 分 成 B=4096 个 8x8 子 块 ;J; 是 第 i 块 的 平均 能 量 , 故 8x8 9:dct2(T);// 将 原 像素 矩阵 进行 DCT 变换 
子 块 的 平滑 度 为 16:compute Fcontrast (i,j) ;// 按 照 2.2 节 公式 
ei=6lx 本 +8x 杞 (15) ”2 ”结合 纹理 复杂 度 和 JND 模型 的 图 像 水 印 算法 
然后 根据 以 下 规则 来 划分 平滑 和 非 平滑 区 : 2.1 水 印 钥 入 过 程 
True， 平滑 区 先 将 水 印 图 像 像素 依次 由 十 进 制 转换 为 二 进 制 。 比 如 : 
Ee 非 平滑 区 ‘9 0 一 00000000，255 一 11111111， 若 水 印 图 像 的 大 小 为 MixN1， 
图 像 中 的 纹理 指 的 是 图 像 中 具有 一 定 随 机 性 、 重 复 性 的 部 则 转换 成 二 进 制 后 的 编码 长 度 为 MixNix8， 将 其 保存 在 密 匙 数 


分 ， 边 缘 指 的 是 图 像 中 具有 -一定 的 对 比 度 和 结构 化 的 部 分 ， 图 ”组 Key[3]-[&.MIINH 中 。 水 印 具体 嵌入 过 程 如 下 

像 的 边缘 区 域 对 比 度 一 般 较 大 、 像 素 之 间 的 相关 性 较 小 ， 因 此 a) 若 原始 图 像 大 小 为 M2xN2， 将 其 进行 32x32 分 块 ， 则 共 
通过 对 比 度 和 性 将 非 平滑 区 划分 纹理 区 和 边缘 区 。 

和 有 mn= ?2x22 图 像 块 ， 对 每 块 图 像 块根 据 式 (1) ~ (4) 分 别 


首先 计算 整 幅 图 像 的 纹理 度 : 32x32 
计算 各 自 的 JHD,COV， 并 根据 式 (9) 求 得 每 块 图 像 块 的 纹理 
复杂 度 f1,f,.…,fh。 
再 计算 8x8 子 块 的 纹理 度 : b) 对 上 步 求 得 图 像 纹 理 复杂 度 按 降序 排序 ， 优 先 从 图 像 复 
高 的 子 块 进行 水 印 嵌入 。 

0) 将 每 个 32x32 子 块 分 成 16 个 8x8 小 块 ， 


56=0, x|D|+0, xCOV (17) 


崔 
洱 


0.=0, x|D,|+0, xCOV, (18) i 
对 每 一 小 块 计 
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录用 稿 


算 JND 值 。 过 程 如 下 : 
QO 计算 每 个 8x8 小 块 JHD,COV， 将 图 像 进行 平滑 、 边 缘 


和 


流 


嵌入 ， 即 每 个 32x32 
32x32 的 水 印 
32x32 图 


纹理 区 域 的 分 类 。 


对 原始 图 像 像素 矩阵 进行 DCT 变换 ; 


名 结合 第 中 步 得 到 的 分 区 结果 , 根据 
像 块 的 JND 值 。 


按照 图 像 纹理 复杂 度 降序 选择 图 像 子 块 进行 水 印 二 进 第 


码 的 嵌入 ， 和 骨 入 规则 如 下 : 


JND(i,j)=0, 不 嵌入 水 印 
i 当前 水 印 编码 为 1: 了 
JND(i, j) 0, 
当前 水 印 编码 为 0: 了 本 


式 (10) 计算 每 个 8x8 


Le 


D(i,j) +dct(i,j) 
D(i,j) -dct(i,)) 


根据 此 规则 ， 每 个 8x8 的 图 像 块 中 5 


的 图 像 块 中 可 嵌入 


图 像 共有 


7 
4 有 >.i= 28 个 位 置 可 
i=l 


6:F=idct2 (DCT) // 将 嵌入 水 印 后 的 DCT 和 珑 阵 进行 逆 变 换 


7:Return F 
2.2 水 印 提取 过 程 
a) 对 原始 


DCT1(i,j), DCT20,)); 


入 位 置 数 


值 , 根据 需要 的 嵌入 位 置 数 


图 像 和 含水 印 图 


像 分 别 进 


b) 根 据 秘 钥 Key[3] 能 够 获得 水 印 图 像 的 大 小 以 及 需要 的 构 
目 ， 即 8xM1xN1, 计算 原始 


图 像 的 纹理 复杂 度 和 JND 
目 、32x32 子 块 复杂 度 排 序 、8x8 子 


行 DCT 变换 得 到 


块 内 允许 的 嵌入 位 置 ( 即 IND#0 的 位 置 ) 获得 每 一 个 水 印 编码 


巾 入 的 位 置 ; 


0) 得 到 位 置 后 ， 根 据 如 下 规则 提取 1 
+JVD 亿 让 , 则 提取 水 印 编码 为 1 
一 JND(i,j), 则 提取 水 印 编码 为 0 
将 其 转换 为 十 进 制 ， 按 


DCT,(i,j)- DCT (i,)) -| 


提取 了 水 印 图 像 的 二 进 制 编码 ， 昨 


28x16=448 个 ， 一 幅 


32x32x8 个 二 进 制 码 , 那么 就 需要 19 块 


的 规定 位 置 按 序 嵌入 水 印 的 二 进 制 编码 。 


经 
算法 2 所 示 。 
结合 纹理 复杂 度 和 JND 模型 的 图 像 水印 柑 入 算法 


d) 对 
过 数据 类 型 的 转换 ， 得 到 含水 印 


图 像 。 


算法 2 
Input: 1)I // 原 始 图 像 ; 

2)SY// 水 印 图 像 。 
Output:F // 含 有 水 印 图 像 的 像素 矩阵 。 


像 块 进行 水 印 的 戏 入 。 所 以 只 要 在 图 像 复杂 度 前 19 块 


联 入 水 印 的 二 进 制 编码 后 的 图 像 算 阵 进行 IDCT 变换 ， 


图 像 水 印 藤 入 算法 如 


1:T1=GLCM(I);// 


2:compute 万 // 按 


原 矩阵 得 到 4 个 方向 的 灰 度 共生 矩阵 
轨 2.1 计算 纹理 复杂 度 排序 


3:compute JND (i,j) ;按照 算法 1 计算 JND 值 


4:SY1[n]=dec2bin(SY);// 将 水 印 各 像素 值 按 序 转 换 成 二 进 制 ， 


一 维 数组 SY1 [n] 


存 入 


5:for num=1,2,...,n // 依 次 租 入 n 个 水 印 编码 


for /,/,…, fis // 按 纹理 复杂 度 依次 散 入 


Cut into temp[8][8];// 按 序 将 各 块 分 割 成 8*8 子 块 


For T=; 250s38 


for j=1,2,...,8 // 对 子 块 各 位 置 进 


if JND (i,j) =0 


j++; // 不 嵌入 
else if SY1i[num]=1; 
DCT(i j)=DCT(i, j)+IND Gi,j) 
else 
DCT(i,j)=DCT(i, j)-IND (i,j) 
end if 
end for // 下 一 列 
end for// 下 一 行 
end for// 下 一 子 块 
end for // 全 部 水 印 编码 嵌入 完毕 


行 操作 


腾 入 的 


具体 信息 : 


序 放 入 MxN 的 矩阵 即 可 得 到 水 印 图 像 。 水 印 提取 的 过 程 如 图 


1 所 示 。 
密 名 
Key[3]= | 瑶 体 图 人 含水 日 像 
[8，M，N] 加 | = 
获取 水 印 图 计算 图 像 纹 得 行 DCT 变换 得 | 进行 DCT 变 换 得 
像 的 大 小 理 复杂 度 到 DCT2 (i,j) ， 到 DCT1 (i, j) 
获取 水 印 嵌入 | 获取 水 印 的 | 
| 的 位 置 | | 
根据 提取 水 印 规 
过 和 二 提 雪 和 到 水 外 国 多 


图 1 水 印 提取 的 过 程 


3 ”实验 结果 与 分 析 


实验 平台 : 本 实验 所 | 
处 理 器 为 酷 害 i5-7200U， 主 频 为 2.6 GHz， 内 存 为 4GB。 仿真 


的 操作 系统 为 


Windows7 旗舰 版 ， 


软件 为 MATLAB 和 codeblock。 实 验 选 择 Lena Barbara、 Lake、 


Boat 四 幅 512x512 标准 
骨 入 水 印 图 像 对 比 


3.1 
4(a) 为 原始 图 


灰 度 图 像 〈.bmp )。 


片 ，(b) 和 (c) 分 别 为 租 入 8 192 bit、32 768 
bit 水 印信 息 后 后 的 图 片 。 


(9) 原 图 


图 2 


(b) 含 8 192 bit 水 印 


嵌入 水 印 前 后 图 像 


(@) 含 3 2768 bit 水 印 


可 见 ， 在 租 入 容量 为 8 192 bit 与 32 768 bit 水 印信 息 后 ， 
像 的 质量 在 视觉 上 几乎 没有 变化 ， 水 印 隐藏 性 较 好 。 四 幅 图 
像 实 验 的 具体 参数 性 能 见 表 2。 


pa 
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表 2 本 文 算法 与 文献 [6JPSNR 的 比较 
PSNR/dB 名 
嵌入 水 印 容量 /bit ”载体 图 像 
文献 [6] 本 文 忒 
fi 
Lena 61.117 64.0496 
Barbara 60.823 62.8909 
8192 
Lake 60.982 63.9762 
Boat 61.736 65.6315 0.5 1 1.5 2 2.5 3 3.5 
志和 实时 4p: 
Average 61.4145 64.1370(+4.43%) 嵌入 容量 /10 bit 
Lena 59.671 61.7252 图 4 Barbara 实验 结果 对 比 
Barbara 55.721 58.0056 
16384 
Lake 58.377 60.9047 Ee 
Boat 60.541 63.3006 64. a 
Average 58.5775 60.9840(+4.11%) 
Lena $4.4561 四 
Barbara ”超过 幅 入 49.7129 Esg. 
32768 了 
Lake 容量 上 限 52.1042 O56- ©. 
Boat 57.5262 54- 
Average 53.4498 52 - © 
505 1 15 2 25 3 35 
要 Re 嵌入 容量 /104bit 
从 表 2 可 见 ， 在 租 入 8 192 bit 水 印信 息 后 ， 本 文 方法 的 平 
图 5 Lake 实验 结果 对 比 
均 PSNR 为 64.137 0 dB， 比 文献 [6] 提 高 了 4.43%; 在 仍 入 16 
384 bit 水 印信 息 后 ， 本 文 方法 的 平均 PSNR 为 60.9840dB， 提 
高 了 4.11%; 而 本 文 方法 在 葵 入 32 768 bit 水 印信 息 后 ， 平 均 68- 
Q Pk % 
PSNR 为 53.449 8 dB， 减 小 了 水 印 嵌 入 带 来 的 图 像 失真 。 Bol 2 
为 了 更 直观 地 说 明 本 文 方法 将 纹理 特性 和 DCT 域 JND 值 ee 
GQ 
联合 鞭 入 水 印 算法 的 优越 性 ， 特 用 折线 图 将 本 文 方法 和 文献 [6] 本 64 芭 
人, ee Ke “i 
单独 利用 纹理 特性 的 方法 着 重 作 比较 ， 四 幅 图 像 各 水 印 容量 骨 ES 62 Se 人 
入 实验 结果 如 图 3~6 所 示 。 @ ~ 
60 - 
58 - 网 
68- SEE 汪汪 © 
S. …@… 本 文 方法 56- 有 
66- 、 一 一 文献 [6] $5 1 1.5 2 2.5 3 3.5 
嵌入 容量 /104bit 
画 图 6 Boat 实验 结果 对 比 
BG 
多 3.2 提取 的 水 印 对 比 
D_ 区 三 ee 
7 展示 了 文献 [6] 和 本 文 方法 提取 的 水 印 (8 192 bit)。 
上 重 
人 wt 二 Wl; 
“© 本 一 
上 r ls r r : ” | | 
S65 1 1.5 2 2.5 3 3.5 月 吧 目 已 HE 
嵌入 容量 /104bit 


图 3 ”Lena 实验 结果 对 上 


t 


原水 印 


提取 的 水 印 
妈 7 ”提取 前 后 的 水 印 


图 


文献 [6] 提 取 的 水 印 
像 
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亲朋 有 
如 图 7 所 示 ， 本 文 算法 提取 出 的 水 印 较 文献 [6] 清 晰 。 文 献 。 [2] 黄 祥 , 杨 式 年 . 度 和 基于 动态 窗口 的 纹理 特征 的 这 感 影像 分 类 
[6] 提 取出 水 印 的 相似 度 为 91.79%， 水 印 的 PSNR 为 37.709 7 [J]. 测绘 科学 技术 学 报 , 2015, 32 (3): 277-281. 
dB, 本 文 方法 提取 的 水 印 相似 度 为 97.55%， 水 印 的 PSNR 为 [3] Haralick R M, Shanmugam K, Dinstein I. Textural features for image 
49.285 1 dB， 从 实验 结果 来 看 ， 在 小 容量 水 印 的 嵌入 提取 中 ， classification [J]. IEEE Trans on Systems Man & Cybernetics, 1973, SMC3 
本 文 的 改进 算法 更 能 准确 地 提取 出 水 印 的 正确 信息 。 (6): 610-621. 
文献 [6] 的 柑 入 规则 限制 了 水 印 的 髋 入 容量 , 该 法 的 散 入 上 [4] 张 浩 , 杨 永 健 , 韩 鸿 葛 . 一 种 基于 DCT 的 序列 分 组 水 印 算法 [J]. 计算 
限 为 16 384 bit， 而 本 文 方法 在 仍 入 32 768 bit 水 印信 息 时 ， 本 机 应 用 研究 , 2014, 31 (11): 3473-3476. 


文 方法 提取 的 水 印 相 似 度 为 87.21%， 水 印 的 PSNR 为 15.7470 [5] Shen X, Shi Z, Chen H. Splicing image forgery detection using textural 


dB， 实验 结 果 见 图 8。 features based on the grey level co-occurrence matrices [J]. Iet Image 


Processing, 2017, 11 (1): 44-53. 


龙 的 Tr 片 1 [6] 李 淑 芝 ， 胡琴， 邓 小 鸿 . 灰 度 共生 引 阵 纹理 特征 选 块 的 可 涝 图 像 水 印 
中 . 光电子。 激光, 2017, 28 (4): 411-418. 
人 人 上 中 (4) 


[7] Nguyen PB, Beghdadi A, Luong M. Perceptual watermarking using a new 


just-noticeable-difference model [J]. Signal Processing Image 


原水 印 提取 的 水 印 
图 8 容量 为 32 768 bit 提取 前 后 的 水 印 图 像 


Communication, 2013, 28 (10): 1506-1525. 
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